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ABSTRACT 


We report on the development of an analytical optical model for scat- 
tering of light due to lightning by clouds of different geometry. We use the 
self-consistent approach and the equivalent medium concept of Twersky to 
treat the case corresponding to outside illumination. Thus, the resulting 
multiple scattering problem is transformed with the knowledge of the bulk 
parameters, into scattering by a single obstacle in isolation. Based on the 
size parameter of a typical water droplet as compared to the incident wave 
length, the problem for the single scatterer equivalent to the distribution 
of cloud particles can be solved either by Mie or Rayleigh scattering theory. 
The super computing code of Wiscombe can be used immediately to produce 
results that can be compared to the Monte Carlo computer simulation for 
outside incidence. 

A fairly reasonable inverse approach using the solution of the outside 
illumination case has been proposed to model analytically the situation 
for point sources located inside the thick optical cloud. Its mathematical 
details are still being investigated. When finished, it will provide scientists 
an enhanced capability to study more realistic clouds. 

For testing purposes, the direct approach to the inside illumination of 
clouds by lightning is under consideration. Presently, we are on the verge 
of obtaining an analytical solution for the cubic cloud. For cylindrical 
or spherical clouds, we need preliminary results of scattering by bounded 
obstacles above or below a penetrable surface interface. 
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I. INTRODUCTION 


We consider the scattering of light (visible or infrared) due to lightning 
by cubic, cylindrical, and spherical clouds. In this report, a typical cloud is 
represented by a statistically homogeneous ensemble of configurations of N 
identical and aligned spherical water droplets whose centers are uniformly 
distributed in its volume V. The radius of a droplet varies from 5 to 15 
microns. The incident light is from outside the penetrable cloud. 

The optical effects of clouds on the light produced by lightning have 
received great interest for many years. Different techniques have been used 
in trying to explain the complicated nature of these effects. In particular, we 
mention the Monte Carlo method which is a computer simulated technique. 
In a Monte Carlo program, we follow the path of the photons emitted into 
the cloud by lightning. A photon is said to be scattered if it escapes from the 
cloud after colliding with the spherical droplets. Otherwise, it is considered 
as being absorbed by the cloud [1 to 8]. 

The two main limitations of the Monte Carlo method [9] are the limited 
number of photons which can be traced through an optically thick cloud with 
finite computer time, and the fact that a Monte Carlo program cannot be 
formally proven to be mathematically correct for a finite number of photons. 
Also, it is very difficult to obtain reliable statistics with the Monte Carlo 
program. 


In this paper, we extend to cloud physics the work done by Twersky [10 to 
12] for single and multiple scattering of electromagnetic waves. We solve the 
interior problem seperately to obtain the bulk parameters for the scatterer 
equivalent to the ensemble of spherical droplets. With the interior solution or 
the equivalent medium approach, the multiple scattering problem is reduced 
to that of a single scatterer in isolation. Hence, the computing methods 
of Wiscombe [13] specialized to Mie scattering can be used to generate 
numerical results to compare with those given by different literatures. 
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II. OBJECTIVES 


The primary objective of this project is to obtain an analytical optical 
model for the scattering of light due to lightning from different point sources 
located inside the cloud. We seek acceptable results in short computer time 
which will take into account the cloud physics and its boundary effects. To 
accomplish this objective, it was necessary to solve the multiple scattering 
problem with outside incidence. By varying the outside angle of incidence, 
we will obtain for a fixed point inside the cloud an orthogonal set of radiative 
solutions. From the same point source inside the cloud, we impose a 
uniformly radiative answer and expand it in terms of the orthogonal set 
of radiative solutions. The expansion coefficients will determine the outside 
incidence corresponding to the uniformly radiative imposed inside solution 
(see figure #1). 


Geometry For Points Sources Inside The Cloud 


(i). The Cubic Cloud 


(ii). The Cylindrical Cloud 




(iii). The Spherical Cloud 
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III. SOLUTION INSIDE THE CLOUD FOR OUTSIDE INCIDENCE 


The solution inside the cloud for outside illumination corresponds to 
the multiple scattering of a plane electromagnetic wave by an ensemble of 
configurations of N identical and aligned spherical water droplets. To obtain 
the solution inside the cloud, we first consider (see figure #2): 

(i). The Single Scatterer In Isolation 



For an incidence plane electromagnetic wave (f> = ke lKx ' r = kt]\ 
and rf being the complex relative index of refraction for the host medium 
inside the cloud, the total outside solution 

xfi=<£+u 0 ( 1 ) 

satisfied the following differential equation obtained from Maxwell’s equa- 
tions after suppressing the harmonic time dependence factor e~~ lu;t 


V x V x ^ = 0, V • ifc = 0. 


(2) 


The solution inside the single spherical water droplet in isolation ip r 
satisfied 


V x V x +kI 


V ; tn — 0, V * Xpin — 0. 


Here, 


n » / // 
k 2 = k y t) — krj rj , 


(3) 

(4) 


with if being the complex relative index of refraction for the medium inside 
the spherical water droplet. 
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Similar to Twersky [11], we have from (1) 


xf = ae‘*‘ r + j/* («i |r - r'|) ,u c (r')| 


with 


u 0 (r') = |h,u Q } 

= j (j 1 * nj ■ (v x u 0 ) - x hj • (n x u Q )j dS (r') 


Here, 


h= 1 + 




1 ). 


and 


h(x) = —, 

IX 

and I being the identity dyadic. 

Asymptotically, for K,r >> l, we can write 

u 0 (r) = h(K 1 r)g(r,'k 1 :a),r-g = o 

and 

g(r)= !< -g(r), 
i, = (i - rf) . 

The scattering amplitude 

g(?,K, : a) = |l 1 e _, ' ?i r ',u 0 (r')J 


(5) 


( 6 ) 


( 7 ) 


(8) 


can also be evaluated from Mie scattering theory. 

The spectral representation of the scattered wave for a single spherical 
water droplet is 


u 0 (r) = — / e'* lc r g (r) dQ (6 c ,ip c ) 

27T J 

c 

with 

r > (r • r') , k 1c = k 1 c t c ( d c , y c ) . 


(9) 
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From general reciprocity relations, we have Twersky [11] 


-9? [a g : a)] = -- - - s 

°0 

with 

4tt (7 0 1 r r | 

= 2 i ’ 6 '*) 

and 

= M |g(r,/c : a)| 2 . 


( 10 ) 


Here, a a and a s are the absorption and scattering cross sections. The 
operator M is the mean over all directions of observation. 

Equations (1) to (10) solve the single scattering problem. 

(ii). A Fixed Configuration N Identical Scatteres (see figure #3) 



We consider a fixed configuration of N identical scatterers with centers 
located by r m(m=1 )2>3 
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The total outside field (r) for an incoming wave <j> is 

N 

*(r) = ,?(r)+£u m (r-r m ) 

m= 1 

U w (r - r m ) = U m (r - r m ; r,, r 2 , r N ), 


( 11 ) 


with 

U m (r — r m ) ~ h (k x |r — r TO |) G m (12) 

for |r — r m | — ► oo. The multiple configurational scattering amplitude is given 
by r i 

G m (f)= {l t e-^ r ',U m }. (13) 

Equivalently, for the scatterer located at t, we use the self-consistent ap- 
proach of Twersky [10] to obtain the total outside configurational field 


Vt (r) = (£(r) + E m u m (r - r m ) 
+U t (r - r<) , 

with 


El = E- 


m^t 


Using (14) and the general reciprocity relation 



(14) 


(15) 


for any arbitrary direction of incidence, we derive as in [11] the self-consistent 
integral equation for the multiple configurational scattering amplitude 


G t (f ) = g t (r, «,) ■ ae’* r ‘ + El / ^ • G "> (?c)e^ c R, “, 


with 

Rtm ~ ** m j ^ 77 d 


1-hS 


( 16 ) 


dil,. 


(iii). Ensemble of Configurations 
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We take the average of (16) over a statistically homogeneous ensemble 
of configurations of N identical and aligned scatterers whose centers are uni- 
formly distributed in the volume of the thick cloud. Using the equivalent 
medium approach and Green’s theorems, we obtain [11] the dispersion rela- 
tion determining the coherent parameters 



+p I [/(R)- i]e- ,S R Ud(R). 

Voc-V 


(17) 


Here, Q is the equivalent scattering amplitude and U is the radiative function 
defined !>y 


|A|R| + U = 0, 
u = J g(f ,r c )-«? (* 1C |K) e , * lcR , 


N 

V 


and the bulk propagation parameter is 


(18) 


K = KjT] 


(19) 


with rj being the bulk index of refraction. 

In equation (17), / (R) — lj is the total correlation function, and 

(boo — (') represents the depleted volume ( the volume of all space less that 
of the exclusion region ). The operator 


{[/*$]} = J [fdng - gd n f]dS (20) 

s 

is the Green exclusion surface operator with outward unit normal from v. 

Equation (17) solves formally the interior problem for the cloud. To 
obtain numerical results, one can apply stationary phase method [12] on 
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(18) and reduce (17) to 


£ = 


K - 


r^j 




and 



J [/(R)-i]e*^- K ) R d(R) 

0 


( 21 ) 


To obtain equation (21), we neglect back scattering in the forward direction 
since the spherical droplets are large tenuous scatterers and the magnitude 
of the single scattering amplitude is evaluated in the forward direction. 
Equation (21) can be solved numerically, iteratively or directly subject to 
explicit restriction on 


g = eg(K,K). 


( 22 ) 


It is important to stress that equation (21) can be subjected to different 
approximations procedures [14] depending on the cloud composition and the 
nature of its boundary. In particular, we mention the asymptotic leading 
term approximation, the generalized Rayleigh’s approximation, and the two- 
space scatterer formalism leading term approximation [11]. 

The two-space scatterer formalism can be used directly in cloud with a 
slab geometry, and it will preserve the self-consistency of the solution. The 
scatterers in the slab geometry will be excited by a coherent wave traveling 
in the equivalent medium but will radiate into free space. 
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IV. THE EQUIVALENT MEDIUM APPROACH 


Now, we can use the results of part (III) to solve the outside illumination 
for clouds with different geometry since the net behavior of the interior 
is known. The multiple scattering problem becomes the scattering by a 
single obstacle (see figure # 4 ). To produce numerical results, we can apply 
Wiscombe’s advanced super computer code specialized to Mie scattering 
Geometry of The Equivalent Medium Approach 

(i). The Cubic Cloud (front view) 



(ii). The Cylindrical Cloud 



(iii). The Spherical Cloud 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The analytical solution for the outside illumination of clouds by lightning 
using the equivalent medium approach has been completed. The multiple 
scattering problem has been reduced to that of a single scatterer in isolation. 
Depending on the size parameter of the cloud particles as compared to the 
wave length of the incident light, either Rayleigh or Mie scattering technique 
can be used to determine Qext, Qscat, and Qbacs. 

With the bulk parameters, we can use Wiscombe’s computer code to ob- 
tain in short computer time, acceptable numerical results for a medium with 
a complex relative index of refraction which is an improvement of Bohren 
[15]. The equivalent medium approach gives naturally the polarizations and 
the angular distributions of photons which escape the cloud surface. 

Due to the complexity of the problem and the time constraint, the 
mathematical details for points sources inside the cloud are still being worked 
out. They will be the subject of more studies in the future. To the best of 
my knowledge, no analytical solution of the problem has yet been published. 
Even partial results will play a major role in describing the transport of 
radiation produced by lightning from point sources located inside the cloud. 

To test the analytical solution, we intend to solve the cubic case directly 
since we know the results of scattering by a penetrable plane interface. For 
cylindrical and spherical clouds however, one has to establish preliminary 
results for the scattering of light by an obstacle above or below a penetrable 
surface interface. 

With the completion of the analytical model, we will have the ability 
and the capability to study the transfer of radiation through realistic clouds. 
That is to say, we will be able to consider clouds made up of water and ice 
particles, of inhomogeneous distribution of diversified constituents; we will 
be able to study the fluctuational changes of polarization of light through 
the cloud due to sudden geometrical deformations of its boundary and the 
frequent realignment of its particles (e.g, needle shape ice particles aligned 
in a strong electric field). 
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